Spin relaxation in n-type GaAs quantum wells with transient spin grating 



M. Q. Wen£] and M. W. WrQ 

Hefei National Laboratory for Physical Sciences at Microscale, 
University of Science and Technology of China, Hefei, Anhui, 230026, China and 
Department of Physics, University of Science and Technology of China, Hefei, Anhui, 230026, Chin^ 

H. L. Cui 

Department of Physics and Engineering Physics, 
Stevens Institute of Technology, Hoboken, NJ 07030, USA 
(Dated: February 1, 2008) 

By solving the kinetic spin Bloch equations, we study the time evolution of the transient spin 
grating, whose spin polarization varies periodically in real space, confined in (001) GaAs quantum 
wells. With this study we can investigate the properties of both the spin transport and the spin 
relaxation at the same time. The Fourier component of the spin signal decays double exponentially 
with two decay rates l/r+ and l/r_. In high temperature regime, the average of these two rates 
varies with the grating wave- vector q quadratically, i.e., (1/7+ + l/r_)/2 = D s q 2 + 1/t s , with D s 
and f s representing the spin diffusion coefficient and the average of the out-of-plane and the in-plane 
spin relaxation times respectively. t± calculated from our theory are in good agreement with the 
experimental data by Weber et al. [Phys. Rev. Lett. 98,076604 (2007)]. By comparing D s with and 
without the electron-electron Coulomb scattering, we calculate the contribution of Coulomb drag to 
the spin diffusion coefficient. With the transient spin grating result, we further reveal the relations 
among different characteristic parameters such as spin diffusion coefficient D s , spin relaxation time 
t s , and spin injection length L s . We show that in the presence of the Dresselhaus and/or Rashba 
spin-orbit coupling, the widely used relation L s = V D s t s is generally inaccurate and can even be 
very wrong in some special cases. We present an accurate way to extract the steady-state transport 
characteristic parameters from the transient spin grating signals. 

PACS numbers: 72.25.-b, 72.25.Rb, 72.25.Dc, 75.40.Gb, 71.10.-w 



I. INTRODUCTION 

Recently a lot of efforts have been devoted to the 
study of spin dynamics and spin transport in semicon- 
ductor nano-structures in order to realize the spintronic 
device* 1 -' 2 -' 3 -'^ In n-type zinc-blende semiconductors, elec- 
tron spins are randomized by the Dresselhaus and/or the 
Rashba spin-orbit coupling (SOC) which acts as an effec- 
tive magnetic field h(k) with its direction and magnitude 
depending on the electron momentum k. 5,6,7 In spatial 
homogeneous system the spin evolution is characterized 
by the spin relaxation time t s which describes the decay 
rate of spin polarization; The steady-state spin transport 
is characterized mainly by the spin injection/diffusion 
length L s j&^ i 10 : 11 while the transient spin propagation 
is characterized by spin diffusion coefficient D s i 12 ' 13 The 
relations among these three parameters and other pa- 
rameters such as momentum relaxation time r p , charge 
diffusion coefficient D c and mobility /i have been actively 
discussed. 

It is understood that spin relaxation/dephasing is in- 
duced by the inhomogeneous broadening due to the SOC 
together with the (spin conserving) scattering ] 14 : 15 The 
scattering provides a channel to speed up the spin de- 
phasing but also slows down the spin dephasing by weak- 
ening the inhomogeneous broadening^ The competing 
effects of scattering have different results on spin de- 
phasing in different conditions, fn weak SOC or strong 
scattering regime, the out-of-plane spin relaxation time 



is expected to be l/r s cx (h 2 (k)r*). ft should be 
noted that r* not only includes the contribution from 
conventional momentum scattering mechanisms such as 
electron-impurity and electron-phonon scattering, but 
also includes the contribution from electron-electron 
Coulomb scattering which does not directly affect the 
charge transport properties. 14 ' 16 : 17 As for the steady- 
state spin injection problem, by assuming that the spin 
dynamics can be separated into two independent pro- 
cesses, spin diffusion and spin relaxation, the spin polar- 
ization is expected to decay exponentially along the in- 
jection direction with decay "rate" , i.e., the spin injection 
length, L s = y/D s T s , where D s and r s are two phenome- 
nal parameters whose relations with other properties are 
yet to be determined j^i^ i 1 1 1 1 However in the presence 
of the SOC, it has been proved that this assumption is 
oversimplified. By solving the kinetic spin Bloch equa- 
tions, it is shown that the effective magnetic field due 
to the SOC alone causes the electron spin to process in 
real space even in diffusive regimei 18 : 19 The spin polar- 
ization varies in the space as e~ x l Ls cos(x/Lq) instead 
of simple exponential decay. The spin injection length 
L s and the spatial oscillation "period" Lq are obtained 
by solving the kinetic spin Bloch equations, 2 — the spin 
transport equations which include the contribution of the 
SOC2i or the linear response theory^ 2 - In the diffusive 
regime, L s and L$ can be expressed by SOC strength and 
momentum scattering timei 20 ' 21 ' 22 As for the spin diffu- 
sion coefficient, it was widely assumed to be the same as 
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the charge diffusion coefficient. Later it was pointed out 
that spin Coulomb drag (SCD), caused by the electron- 
electron Coulomb scattering, should suppress the relative 
motion of electrons with different spins and thus reduce 
D 13 - 23 - 24 - 25 

A direct measurement of spin diffusion coefficient 
can be carried out by transient spin grating (TSG) 
experiments! 12 ' 13 i 26 ' 27 With the assumption that spin dif- 
fusion and spin relaxation are independent of each other, 
the decay rate of TSG was written as T q = D s q 2 + l/r s , 
where q is the wave- vector of the spin grating, 12113 How- 
ever, from the lesson of steady-state spin injection one 
learns that the spin diffusion and relaxation are not sep- 
arable even in the diffusive regime. This is further justi- 
fied by the fact that the decay of TSG can be fitted to 
a double-exponential form instead of single exponential 
one.— It is therefore still a question of how to get spin 
diffusion coefficient through TSG experiments. 

In this article we study the temporal evolution of the 
TSG by solving the kinetic spin Bloch equations. This 
paper is organized as following: In Sec. [IT] we first set up 
the kinetic spin Bloch equations and apply them to solve 



the TSG problem. Then we show the analytical solution 
of simplified equations and the numerical results of the 
full kinetic spin Bloch equations. In Sec. lIIII we study the 
relations among the spin relaxation time r s of a spatially 
homogeneous system, the spin injection length L s and 
the spatial oscillation period Lq of the steady-state spin 
transport as well as D s of the transient spin transport. 
We conclude in Sec. IV. 



II. KINETIC SPIN BLOCH EQUATIONS 

The system we study is the electron gas confined in 
a (001) GaAs quantum well (QW) with width a grown 
along the z-axis. We assume that the well width is nar- 
row enough so that only the lowest subband is occu- 
pied. With the help of the noncquilibrium Green func- 
tion method^ one can write down the kinetic spin Bloch 
equations ^ 18 ' 19 by using gradient expansion and gener- 
alized Kandanoff-Baym ansatz: 



dpk(M) , 
dt 

_ dp k (x, t) 
dt 



eE(x) 



dpk{x,t) k x dp k (x,t) 



dk x 



Ox 



iKffMsB + h(k)) • a/2 + £„ F (a;, k, t),p k (x, t)] 



(1) 



Here we assume that the transport direction is along the 
x-axis. p(x, k, t) is the density matrix whose diagonal el- 
ements /ko-(^,i) represent the electron distribution func- 
tions with spin a{— ±1/2) and momentum k = (k x ,k y ) 
at position x. The off-diagonal elements stand for the 
spin correlations between spin- up and -down electrons. 
The second and third terms of Eq. (fTJ) correspond to 
the drift driven by the electric field E(x), determined 



by the Poisson equation, and the diffusion of electrons, 
respectively. The fourth term describes the spin preces- 
sion around the total magnetic field which is composed 
of the external magnetic field B, the effective magnetic 
field h(k) due to the SOC as well as the one from the 
Hartree-Fock term of the electron-electron Coulomb in- 
teraction £uf(x,\s., t). h(k) contains the Dresselhaus and 
the Rashba terms££>£ 



h(k) = f3(-k x cos29 + k y sm29,k x sm29 + k y cos29,0) 
k 2 ~ k 2 

+7( X V sin2 6> + k x k y cos 29)(k y , -k x ,0) + a(k y , -k x , 0) , 



(2) 



where 9 is the angle between x-axis (the spin injec- 
tion/diffusion direction) and the (100) crystal axis^ 
(3 = -fir 2 /a 2 with 7 being the Dresselhaus coefficient^ a 
represents the Rashba parameter which depends on the 
electric field along the growth direction of the QW. The 



scattering term 



dpk{x,t) 



at 



includes all the scattering, i.e., 



the electron-impurity, the electron-phonon and most im- 
portantly the electron-electron Coulomb scattering. It 
is noted that in our calculation the electron-electron in- 
teraction is treated beyond the Hartree-Fock approxima- 
tion. The expressions for the Hartree-Fock and the scat- 
tering terms are given in detail in Ref. [30I 
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The kinetic spin Bloch equations describe the spin dy- 
namics in the presence of drift, diffusion and spin pre- 
cession. By choosing some specified initial and bound- 
ary conditions, one can obtain the evolution of the spin 
signal in time and real space by solving these equa- 
tions for different systems. In order to study the TSG, 
the initial spin polarization of the electrons is chosen 
to be a sinusoidal wave along the cc-direction P{x) = 
(N r (x)-N l (x))/(N r (x) + N L (x)) = P sm(x/L) but uni- 
form along the y-axis, where L = 2n/q is the spatial 
period. Using periodical boundary condition, one only 
needs to study the dynamics in one period of the space 
regime. 

Unless one makes some simplifications, the kinetic spin 
Bloch equations are too complicated to be solved analyti- 
cally. In this paper, we first present an analytical solution 



in the diffusive regime using simplified equations. This 
solution can only provide an intuitive vision of the TSG 
dynamics. We then present the numerical solution of the 
full kinetic spin Bloch equations. 



A. Simplified Solution 

By neglecting the Hartree-Fock term, the inelastic 
scattering such as the electron-phonon and the electron- 
electron Coulomb scatterings and the coupling to the 
Poisson equation, one is able to rewrite Eq. (fTJ) in dif- 
fusive regime when the scattering is strong enough, by 
using a similar method for calculating the spin relaxation 
as in Refs. |5||6|, 



dS(q,t) 



Of 



Dq 2 S(q,t) +iqh x S(q,t) 



/ (a 2 +£j 2 -2a/3sin20)(fc 2 ri) 
' +<(7fc 3 ) 2 T 3 )/16 

2a/3(/c 2 Ti)cos2(9 



2a/3(fc 2 n)cos26i 

(a 2 +/3 2 + 2a/3sin20)(fc 2 Ti) 
+ <(7fc 3 ) 2 T3>/16 
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2(Q 2 +/3 2 )(fc 2 n) , 

+ <(7fc3) 2 T 3 }/8 / 



S(q,t) = 



(3) 



Here S(q, t) — j d 2 kTr{p(g, k, t)cr}d9 is the Fourier com- 
ponent of the spin density function. The second term 
represents the diffusion term with the diffusion constant 
D = (fc 2 ri/2m 2 >. 1/n = cos(W)d6/2ir with 

r(fc, 9) being the momentum relaxation time due to the 
electron-impurity scattering. The third term is the spin 
rotation caused by the non-vanishing net effective mag- 
netic field 

h = (fc 2 Ti/m) (-/3 cos 26,(3 sin 29 - a, 0) (4) 

due to the DP effect and the diffusion. The last 
term is the spin relaxation caused by the DP ef- 
fect and the spin conserving scattering. For a sys- 
tem not far away from the equilibrium, (• • • ) = 
/ • • • df(e k )/de k d 2 k/ J df(e k )/de k d 2 k with f(e) being 
the Fermi distribution function. It is noted that sim- 
ilar results at low temperature have been obtained by 
different approach recently.— 

In the presence of both the Dresselhaus and the Rashba 
terms, the spin relaxation rates are highly anisotropic. In 
the weak SOC regime, the in-plane spin relaxation rates 
are characterized by two decay rates ((a ± $) 2 k 2 Ti)/2 + 
7 2 (fc 6 T3)/32, corresponding to the spin relaxation along 
the characterized directions (110) and (110) axises re- 
spectively. While the out-of-plane relaxation rate is 
the sum of these two rates l/r s = ((a 2 + /3 2 )fc 2 r 1 ) + 
7 2 (fc 6 T3)/16. Note that, for the quasi-two-dimension sys- 
tem, the presence of the cubic Dresselhaus term mod- 



ifies the coefficient of linear Dresselhaus term f3 to be 
(3 — (3 — jk 2 /4. This relation together with Eq. ([3]) are 
briefly derived in the Appendix A. Without the spin ro- 
tation term, the evolution of the TSG is characterized 
by the decay rate Dq 2 + 1/t s . However, the spin rota- 
tion mixes the in-plane and out-of-plane spin dynamics. 
As a result, the evolution is usually described by three 
relaxation rates. There are two special cases where the 
evolution can be described by two relaxation rates. The 
first one is in the system where only one of the Dres- 
selhaus and Rashba terms is important so that the two 
in-plane spin dynamics become identical. The other is 
the spin injection/diffusion along (110) and (110) axises 
where the net effective magnetic field h is parallel to one 
of the characteristic directions of the in-plane spin dy- 
namics. Thus it can only mix the remaining in-plane spin 
with the out-of-plane spin dynamics. The TSG evolution 
of these cases has a double-exponential form, 

S z (q,t) = S z {qMKz- t,T+ + X-e- t/T -) , (5) 
with relaxation rates 

r± = — = Dq 2 + - — + - ± — Wl + f-^ , 

T± 2 T S 1 T s 2t s2 y 'si 

in which 
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Here r s i {t S 2) is the spin relaxation time of the in-plane 
spin which mixes (does not mix) with the out-plane spin 
due to the net effective magnetic field. For spin injec- 
tion/diffusion along (110) axis, r sl = ((a — $) 2 k 2 Ti)/2 + 
7 2 (fc 6 r 3 )/32 and r' sl = ({a - f3) 2 k 2 r 1 )/2. For spin injec- 
tion/diffusion along (110) axis, t s i — ((a + j3) 2 k 2 Ti}/2 + 
7 2 (fc 6 T 3 )/32andr^ = {(a+Pfk 2 ^) /2. In the long wave 
length limit (q < 1), T+ ~ l/r s + (1 + 4T s2 /T' sl )Dq 2 and 
T_ ~ l/r s i + (1 — 4r s2 /Tg 1 )£>q 2 become quadratic func- 
tions of q, roughly correspond to the out-of-plane and the 
in-plane relaxation rates respectively. In general both 
of these two decay rates [Eq. ©] are no longer simple 
quadratic functions of q. If one uses the quadratic fit- 
ting to yield the spin diffusion coefficient, one either gets 
larger (for T + ) or smaller (for T_) values than the true 
spin diffusion coefficient. The accurate way to get the 
information of spin diffusion coefficient should be from 
the average of these two rates 

r - (r+ + r_)/2 = Dq 2 + (l/r s + l/r sl )/2 , (8) 

which differs from the current widely used formula by 
replacing the spin decay rate by the average of the out- 
of-plane and in-plane ones. 
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FIG. 1: (Color online) Spin relaxation times r+ (open) and r_ 
(filled) vs. spin grating wave vector q for three different tem- 
peratures T = 120 (red boxes), 150 (blue triangles) and 295 
K (black circles). The inset is the ratio t_/t+ as a function 
of temperature for q = 0.63 x 10 4 cm -1 . 



B. Numerical Results 

The spin diffusion coefficient obtained by the simpli- 
fied kinetic spin Bloch equations does not include the 
contribution of the spin Coulomb drag since the electron- 
electron Coulomb scattering is neglected. Moreover, the 



simplified equations are derived with only the elastic scat- 
tering. In the equations all of the relaxation times T\ 
that appear in the diffusion coefficient D, net effective 
magnetic field h as well as the spin relaxation matrix 
are the same. For the inelastic electron-phonon scatter- 
ing, it is still possible to write down Eq. ^ with same 
relaxation time t% by the elastic scattering approxima- 
tion. This approximation is valid only for the electron- 
acoustic phonon scattering at high temperature and is 
not valid for the electron-LO phonon scattering. When 
the Coulomb scattering is important, it gives different 
contributions to these relaxation times. On the one hand, 
the relaxation time t\ in the spin relaxation matrix is af- 
fected by the whole Coulomb scattering. On the other 
hand, since the Coulomb scattering among the same spin 
specie does not change the motion of the center of mass, 
it does not directly affect the spin diffusion coefficient. 
Therefore the relaxation time t% in the spin diffusion co- 
efficient is only affected by part of the Coulomb scatter- 
ing. For n in the net effective magnetic field, it is even 
more complicated to analyze the role of the Coulomb 
scattering since it is the joint result of diffusion and spin 
precession. Therefore the relaxation times T\ in the three 
different terms should be different in the present of the 
Coulomb scattering. It would be extremely difficult to 
get the analytical results when the Coulomb scattering 
is taken into account. In order to study the drag effect, 
we numerically solve the full kinetic spin Bloch equations 
Eq. (TT]), with all the scattering explicitly included. The 
numerical scheme is laid out in detail in Appendix B. 
For simplicity, we first consider the spin diffusion along 
(100) axis in a symmetrical QW in which the Rashba 
term vanishes. In our calculation, the SOC strength, the 
electron and impurity densities, and the QW width are 
chosen to be 7 = 11.4 eVA 3 ^ 1 N = 7.8 x 10 11 enr 2 , 
Ni = 1 x 10 11 enr 2 and a — 12 nm respectively. The 
material parameters are listed in detail in Refs. l3Ct 

Our numerical results show that the temporal evolu- 
tion of TSG can not be fitted by a simple exponential 
function with desirable accuracy. However, if we use 
double-exponential function, the accuracy can be im- 
proved more than one order of magnitude. This justi- 
fies that the temporal evolution of the TSG has indeed 
double-exponential form in high temperature regime. In 
Fig. 1 we present the relaxation times of TSG as a 
function of the grating wave-vector q for temperatures 
T = 120, 150 and 295 K. It is seen that r + decreases 
monotonically as q increases while r_ has a peak at some 
small qo. The wave- vector of the peak red-shifts when 
the temperature increases. In the inset of Fig. 1, we 
show the ratio r_/r + as the function of temperature for 
fixed q. One finds that the ratio decreases with temper- 
ature. Our results are in contrast to the predictions of 
earlier theoretical works that the ratio of these two de- 
cay rates and the position of the peak depend only on 
the SOC and material parameters, but do not vary with 
the temperaturei 22 i 32 Recent experiment showed that the 
ratio of these two decay rates are indeed decreases with 
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the increase of the temperature^ It is understood that 
the temperature dependence of the peak position and ra- 
tio between two decay times originate from the cubic k- 
term in the Dressclhaus effective magnetic field. In earlier 
works, it was assumed that only the linear term is im- 
portant. However, in wide quantum wells with high elec- 
tron density, the cubic term becomes important. Once 
the cubic term is considered, the peak qo moves from 
go = vl5m/3/2 (which is independent of temperature) to 
about q Q = vT5m/3/2 = VT5m(/3 - 7 /c 2 /4)/2. Since (fc 2 ) 
increases with temperature, thus qq decreases. The tem- 
perature dependence of r + /r_ also originates from the 
contribution of the cubic Dressclhaus term. 
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FIG. 3: (Color online) F = (r+ + r_)/2 and Ar = (r+ - 
r_)/2 vs. q at T — 295 K. Open boxes/triangles are the 
relaxation rates F + /_ calculated from the full kinetic spin 
Bloch equations. Filled/open circles represent F and Ar re- 
spectively. Noted that the scale for Ar is on the right hand 
side of the frame. The solid curves are the fitting to V and 
AF respectively. The dashed curves are guide to eyes. 



FIG. 2: (Color online) Spin relaxation times t± vs. temper- 
ature for (a) high-mobility sample with q = 0.58 x 10 4 cm" 1 
and (b) low-mobility sample with q — 0.69 x 10 4 cm" 1 . The 
dots are the experiment data from Ref. l27l . 

Due to the natural of the numerical calculation, it is 
not possible to show the results with all possible param- 
eter combinations. However our qualitative conclusions 
are valid not just for this particular parameter set but for 
a large range of parameters. In order to check the quan- 
titative accuracy of our numerical calculation we further 
show the spin relaxation times as functions of temper- 
ature together with the experimental data from Ref. [27J 
in Fig. 2. Curves in Fig. 2(a) are the theoretical spin re- 
laxation times of high mobility (// = 1.5 x 10 4 cm 2 /Vs) 
sample with q — 0.58 x 10 4 cm -1 , while the dots are the 
experimental data from Ref. H?]- Figure 2(b) shows the 
result of low mobility sample (fj, = 3.5 x 10 3 cm 2 /Vs) with 
q = 0.69 x 10 4 cm -1 . In the calculation we use the finite 
square well assumption^ All the parameters we use are 
chosen to be the experimental value if available, eg. the 
grating wave-vector, the electron density, the quantum 
well width, and the impurity concentration determined 
from the mobility. The only adjustable parameters are 
the spin-orbit coupling coefficients 7 and a. In the cal- 
culation, 7 is chosen to be 11.4 eVA 3 and 13.8 eVA 3 for 
the high and low mobility samples respectively and a is 
set to be 0.3/3, close to the choice in Ref. H^- One can see 
from the figure that our theoretical results are in fairly 
good agreement with the experiment data. 

In Fig. 3, we plot the decay rates T± = l/r± and 
their average T = (T + + r_)/2 and difference Ar = 
(r+ - T_)/2 as functions of q at T = 295 K. The de- 



cay rates T± fit poorly with a quadratic function of q. 
In contrast, the average decay rate T fits pretty well by 
the function T = D s q 2 + 1/t' s . The resident error of the 
quadratic fitting for T is two orders of magnitude smaller 
than those of r± . Moreover we find that t' s is very close 
to 4t s /3, inverse of the average of the in-plane and out- 
of-plane spin relaxation rate. For example, at T = 295 K, 
t' s is about 107.8 ps compares to 4r s /3 = 111.6 ps with 
t s calculated by solving the kinetic spin Bloch equations 
for spacial uniform system with the same parameters, 1 - 
Inspired by Eq. |(5J), the coefficient of the quadratic term 
D s can be reasonably assumed to be the spin diffusion 
coefficient. In this way one can calculate the spin diffu- 
sion coefficient with the effect of Coulomb drag included. 
The difference of T + and T_ fits well as a linear function 
of q. The linear coefficient of Ar is about 2^J D s /t' s1 
from Eq. ([6]). 

In Fig. 4 we present the spin diffusion coefficient cal- 
culated in the above mentioned method as a function of 
temperature. For comparison, we also include the charge 
diffusion coefficient, which is calculated by solving the 
kinetic spin Bloch equations with the initial condition 
being the charge gradient instead of the spin gradient. It 
is clearly seen from the figure that D s < D c . For spin- 
unpolarized charge diffusion, the electrons move along 
the same direction and the Coulomb scattering does not 
change the center-of-mass motion, therefore it does not 
change the charge diffusion coefficient directly. However, 
in the spin-polarized transport, spin- up and -down elec- 
trons move against each other and the Coulomb scatter- 
ing therefore slows down the relative motion of these two 
spin species and reduces spin diffusion coefficient. This 
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FIG. 4: (Color online) Diffusion coefficient as a function of 
temperature. Solid circles: Spin diffusion constant D s with 
Coulomb drag; Open circles: Charge diffusion constant D c . 



is the so-called spin Coulomb drag effect i 23 ' 24 ' 25 From 
the figure one can tell that, in the temperature regime 
we study, as the temperature increases, both spin and 
charge diffusion coefficients decrease and their difference 
also decreases. Therefore the Coulomb drag is stronger 
in the low temperature regime. However, even at room 
temperature the Coulomb drag is still strong enough to 
reduce the diffusion coefficient by 30 %. These results 
quantitatively agree with those of Refs. I23H24I 

It is also pointed out that the reduction of spin dif- 
fusion coefficient mostly comes from the Coulomb drag. 
The SOC only has slightly effect on the diffusion coeffi- 
cient since the SOC is very small compared to the Fermi 
energy. The numerical result shows that removing the 
SOC only changes spin diffusion coefficient up to one 
tenth percent for the system we studied. 



III. STEADY-STATE SPIN INJECTION 

In this section we discuss how to obtain the steady- 
state spin injection information from the TSG signal, i.e., 
to find out the relation between the steady-state spin 
injection length L s and spatial spin oscillation "period" 
Lq and the spin diffusion constant D s together with the 
spin relaxation time t s . 

We first show from the simplified solution presented 
in Sec. Ill Al that the steady-state spin injection can be 
extracted from the TSG signal by integrating the TSG 
signal Eq. ^ over the time from to oo and the wave- 
vector from — oo to oo. Form Eqs. ^ and ([7]). the in- 
tegrated TSG reads S z (x) = S z {0)e~ x / L ° cos(x/L + ip) 



L = {JD 2 s t s1 t s I cos- 



(9) 
(10) 



In these equations 

cos = Vm(4/^i - l/r s i - l/r,)/2. (11) 

It is noted that if one only considers the Rashba term or 
the linear Dresselhaus term, t' s1 — t s \, one then recovers 
the result L s = 2\/D s t s from linear response theory^ 
It is seen that the spin precession actually prolongs the 
out-of-plane spin injection length by mixing the fast de- 
cay of the out-of-plane spin with the slow decay of the 
in-plane spin. It is further noted from Eq. ([9]) that the 
spin injection length L s is generally larger than \/D s t s . 
The only exception is when the spin injection is along 
(llO)-direction for the QW with equal linear Dresselhaus 
and Rashba spin-orbit couplings. In this case 6 = 7r/4, 
therefore the net effective magnetic field [Eq. Q] van- 
ishes. Consequently the in-plane and out-of-plane spin 
modes do not mix. 

For most of the cases, L s is in the same order of \/D s t s , 
although the former is usually larger. However, there are 
some special cases where <p — ► 0, L s can be orders of 
magnitude different from y/D s T s . Specifically, according 
to Cheng et al.^ when the spin injection/transport di- 
rection is along (110) in (001) QW, L s becomes larger 
and larger as a approaches (3, regardless of the direction 
of spin polarization. At the limit of a = (3, the spin injec- 
tion length trends to infinity when the cubic Dresselhaus 
term is ignored and the spin oscillates with a spatial pe- 
riod of 27r/(2m/3).— As t s is finite and D s changes little 
for different grating directions, there is no way to obtain 
infinite L s from y/D s T s . 

The infinite injection length and the finite oscillation 
period can also be understood from the TSG point of 
view. Without the cubic Dresselhaus term, l/r S 2 = 
(a — /3) 2 (k 2 Ti) /2 approaches zero and l/r s i — * l/r s = 
2(3 2 {k 2 Ti) when a — > /3. Consequently t± = (yDq ± 
1/v^) -2 wnen ce — p. It is then straightforward to see 
that r_ becomes infinite provided q = qo = 1/ \fDr s = 
2m 0. Therefore the steady-state spin injection along 
(llO) axis is dominated by this non-decay TSG mode 
which is responsible for the infinite spin injection length 
and the spatial oscillation period 2ir/qo = 2ir / (2mP) . 

When the electron-phonon and the electron-electron 
Coulomb scatterings are taken into account, cf> should be 
revised accordingly. Unfortunately, there are no analyti- 
cal expressions for r s , t s \ and t' s1 when all the scatterings 
are included. Nevertheless, one can directly use the TSG 
result to calculate the static injection parameters. The 
numerical result indicates that the two decay rates obey 

r± = D s q 2 + 1/t' s ± (cq + d) (12) 

where c and d are the fitting parameters to 

Ar = cq + d . (13) 
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Q 


Injection 
direction 


D s (cm 2 /s) 


t's (PS) 


c (/im/ps) 


d (ps" 1 ) 


L J S (jmx) 


L (/im) L s Gum) L (fim) 


Q 


= 


(100) 


246 


107.8 


0.02 


3.5 x 10~ 4 


2.27 


2.44 2.28 2.46 


Q 


=0 


(110) 


243 


49.7 


9.4 x 10" 3 


0.012 


0.92 


5.1 1 4.9 


Q 


= 


(110) 


250 


28.1 


0.058 


-2.7 x 10~ 3 


2.7 


0.86 2.1 0.9 



TABLE I: Comparison of spin injection length L s and spatial oscillation length Lo along different injection directions at T = 295 
K from different approaches. The superscripts "S" and "T" stand for L s and Lo obtained from the injection calculation in the 
steady state and those from TSG parameters by using Eqs. (|14[1 and (|15l) respectively. 



Again by integrating the TSG signal over the time and 
wave- vector q with above fitted T±, one obtains 

L s = 2D s /^/\c 2 -4D s (l/T>-d)\ , (14) 
Lo = 2D s /c. (15) 

We stress that these two equations give the right spin 
injection length and the spin oscillation period in the 
presence of the SOC. From the experiment point of view, 
one can monitor the time evolution of TSG with different 
wave- vectors q and obtain the corresponding decay rates 
r± . From the g-dependence of the decay rates, one can 
calculate the spin injection length and spin oscillation 
period from Eqs. flU) and (JTHJ. It is noted that Eq. (fl4|) 
naturally gives the infinite inject ion len gth in the special 
case 2 ^ discussed above as c = 2 y/D s /r s and d = so that 
the denominator in Eq. (fT4)) tends to zero. In contrast, 
\/D s t s always remains finite. 

In order to check that the accuracy of this approach 
[Eqs. (fT4")) and (|15p]. we compare L s and Lq obtained 
from directly numerically solving the kinetic spin Bloch 
equations [Eq. |T])] for spin injection in the steady state 
as described in Refs. |20||29|. with those from the TSG 
approach. In Table [J we list L s and Lq obtained from 
the TSG signal and from the steady-state solution of the 
kinetic spin Bloch equations for spin injection direction 
along (100) axis in symmetrical (001) QW, as well as 
(110) and (110) axes in asymmetrical (001) QW with 
equal Dresselhaus and Rashba coupling. For the sake of 
clarity, we also list the corresponding fitted parameters 
t' s , c and d in the table. It is noted that in the calcula- 
tion the cubic Dresselhaus term is included. One can see 
from the table that even though the injection and the 
oscillation lengths at different conditions are quite dif- 
ferent, the spin diffusion coefficients D s are almost the 
same. Their differences are within the numerical fitting 
error. This confirms that the SOC is too small compared 
to the Fermi energy to affect D s . Moreover, for all the 
cases we study, Lo and L s obtained from these two meth- 
ods agree with each other very well. Although we should 
point out that the accuracy of Lo is higher than L s due 
to the numerical error analysis. Since c 2 — 4D s (1/t' s — d) 
is close to 0, its numerical error is much larger than c. 
Therefore L s calculated from TSG is not as accurate as 
Lo- 



IV. CONCLUSION 

In conclusion, we study the evolution of TSG by solv- 
ing the kinetic spin Bloch equations with periodic bound- 
ary condition. From the simplified equations, we analyti- 
cally show that due to the spin precession the spin signal 
decays double-exponentially instead of simple exponen- 
tially, even in the diffusive regime. Neither of these two 
decay rates are simple quadratic functions of the grat- 
ing wave- vector q. However their average value depends 
quadratically on q and the corresponding coefficient of 
the quadratic term is the right diffusion coefficient D s . 
Therefore it is more accurate to yield the diffusion coeffi- 
cient from q dependence of the average of the two decay 
rates. 

We further show that the corresponding solution from 
the full kinetic spin Bloch equations which include all of 
the scattering mechanisms, especially the Coulomb scat- 
tering, is also in the form of the double exponential decay. 
From the (/-dependence of the average decay rate, one can 
calculate the spin diffusion coefficient with the Coulomb 
drag effect included. It is shown that the Coulomb drag 
effect is stronger in low temperature and decreases with 
the increase of temperature. However, even at room 
temperature the Coulomb drag is still an important fac- 
tor which reduces the spin diffusion coefficient markedly 
compared to the charge diffusion coefficient. 

We also show that by using the TSG result one can 
obtain the characteristic steady-state transport parame- 
ters, such as the injection length and the spatial oscil- 
lation length, from diffusion coefficient D s and the spin 
relaxation times. We point out that in a system with 
the SOC, the wildly adopted relation L s = \/D s t s in the 
literature is generally quantitatively inaccurate and can 
be even qualitatively wrong for some special cases. The 
accurate way to extract the spin injection and spatial 
oscillation lengths directly from the TSG decay rates at 
different grating wave-vector, which are experimentally 
measurable, is proposed. We believe this investigation 
is important in fully understanding the TSG signals in 
experiment. 
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dl} can be derived by expanding Eq. (fT]) in angular 
momentums.-^ Neglecting the Hartree-Fock term, the 
inelastic scattering and the electric field, the Fourier com- 
ponent of the Z-th order of density matrix obeys the fol- 
lowing equation: 



APPENDIX A: DERIVATION OF SIMPLIFIED 
SOLUTION 



By using a similar method for calculating the spin 
relaxation due to the Dresselhaus effect, Eqs. © and 



dpijq, k, t) _ ikq ( 
dl 2m 

pi(q,k,t) 



n 



f p l+1 (q,k,t) +pi-x(q,k,t) j 



LL L — m 

(k) ■ cr,p m (q,k,t) 



(Al) 



where pi (q, k, t) = J e ^ pk(q, t)d<f)/2ii and k = (k cos 0, fcsin0, 0). When both the Dresselhaus and the Rashba terms 
are taken into account, 



h(k) = 7fc((7r/a) 2 — fc 2 /4)(— cos0,sin 0, 0) + 7fc 3 /4(- cos 30, - sin 30, 0) 
+afc(sin 0, — cos 0, 0) , 
= /3k(— cos 0, sin 0, 0) + afc(sin0, — cos 0, 0) + 7fc 3 /4(— cos 30, — sin 30, 0). 

I 



(A2) 
(A3) 



Therefore there are four effective magnetic field compo- 
nents which do not vanish, namely h±i(k) and h±3(k). It 
should be noted that, in quasi-two-dimensional system, 
the ±1 components of the Dresselhaus term are modified 
by the cubic term (3 = ^[(n/a) 2 - fc 2 /4] = [3 — 7& 2 /4. 
When the scattering is strong, one can drop the terms 
with ti higher than the first order and rearrange Eq. (|A1[) 
to obtain Eq. ©. The only difference is that in the spa- 
tial inhomogeneous system, we have additional terms re- 
lying on the wave-vector q. These additional terms give 
rise to the second and third terms in Eq. (3). 



APPENDIX B: NUMERICAL SCHEME 

In order to solve the kinetic spin Bloch equations nu- 
merically, one has to discretize the real space, the mo- 
mentum space as well as the time. The real space is 
divided into segments with equal length. The momen- 
tum space is divided into grids of equal energy and an- 
gular differences ^ The second order up- wind differen- 
tial scheme is applied for the diffusion term and the drift 



term. The former reads 



k x df(x) 
m dx 



k x 3f(x)-4f(x-Ax)+f(x-2Ax) 

Tri 2/S.x 
k x 3f(x)-4:f(x+Ax) + f(x+2Ax) 
m 2Ax 



*«. > o 

m 

fex < Q 



(Bl) 

The boundary condition for x is chosen to be the periodic 
one when we calculate the TSG problem or fixed when 
we calculate the steady-state injection problem.— In the 
energy {£) and angular (0) space, the drift term reads 



eE{x 



dg(k x , k y ) eE(x) ( d 



dk x 



— [2V£ cos 



'2m 
d g(S,9)shx9 



89 



(B2) 



Similar to the diffusion term, one can easily write down 
the second order up-wind differential schemes for £ and 
9 respectively. The numerical schemes for the spin pre- 
cession and the scattering terms are laid out in detail in 
Refs. [USE 

We apply the third order semi-implicit Adams- 
Bashforth scheme for the time differential 3 ^, to achieve 
higher accuracy in temporal evolution. This scheme also 
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saves CPU time. The differential scheme is then given 
by 

p k (x,t + At) -p k (x,t) 5 

^ = ^F[p k (x, t + At)} 

23 

-F[p k (x, t)} + 3/4F[ Pk (x, t - At)} + —G[p k (x, t)} 

16 5 
--G[p k (x,t-At)] + -G[p k (x,t-2At)} . (B3) 

Here F[pi L (x,t)} denotes the drift and diffusion terms 
(the second and third terms) in Eq. ([1]) and G[p k (x,t)} 
stands for the spin precession and the scattering terms 
(the fourth and fifth terms) in Eq. ([1]). The implicit part 
of the equation is solved by Jacobian-free Newton-Krylov 



algorithm^ 

The accuracy of the numerical scheme used in this pa- 
per is higher than the one used in our previous works.— ^ 
The main numerical errors come from the drift and diffu- 
sion parts since the grid size of real space and momentum 
space is limited by the computing power. We find that 
the accuracy of the temporal evolution does not change 
the result of the final steady-state spin transport too 
much. It is therefore expected that the present numer- 
ical scheme and the previous one give very close results 
on the steady-state transport properties. However, the 
present scheme enables us to also study the time sensitive 
phcnomcnons such as the TSG to a sufficient accuracy. 
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